Diffusion of Finite-Sized Hard-Core Interacting Particles In a One-Dimensional Box - 

Tagged Particle Dynamics 
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We solve a non-equilibrium statistical mechanics problem exactly, namely, the single- file dynamics 
of N hard-core interacting particles (the particles cannot pass each other) of size A diffusing in a 
one dimensional system of finite length L with reflecting boundaries at the ends. We obtain an exact 
expression for the conditional probability density function pT(yr ,t\yr,o) that a tagged particle T 
(T = 1, —,N) is at position yr at time t given that it at time t = was at position j/r,o- Using 
a Bethe ansatz we obtain the JV-particle probability density function, and by integrating out the 
coordinates (and averaging over initial positions) of all particles but particle T, we arrive at an 
exact expression for pr{yT, t\yr,o) in terms of Jacobi polynomials or hypergeometric functions. 
Going beyond previous studies, we consider the asymptotic limit of large N, maintaining L finite, 
using a non-standard asymptotic technique. We derive an exact expression for pr (yr , t\yr ,o) for 
a tagged particle located roughly in the middle of the system, from which we find that there are 
three time regimes of interest for finite-sized systems: (A) For times much smaller than the collision 
time t <C t co ii = l/(g 2 D), where g = N/L is the particle concentration and D the diffusion constant 
for each particle, the tagged particle undergoes normal diffusion; (B) for times much larger than 
the collision time t 3> t co u but times smaller than the equilibrium time t <C T cq = L 2 /D we find a 
single-file regime where pr (yr, t\yr,o) is a Gaussian with a mean square displacement scaling as t 1 ^ 2 ; 
(C) For times longer than the equilibrium time t 3> T cq , pr (yr , t\yr ,o) approaches a polynomial- 
type equilibrium probability density function. Notably, only regimes (A) and (B) are found in the 
previously considered infinite systems. 



I. INTRODUCTION 



Recent development of single fluorophore tracking 
techniques allow experimental studies of the motion of 
particles in cellular environments with nanometer reso- 
lution [l[. The cell interior represents a crowded envi- 
ronment, in which the motion of an individual particle 
is strongly effected by the presence of other particles. 
Crowding affects, for instance, the folding of proteins, 
diffusional motion 0, H| as well as rates of biochemical 
reactions 0, [f|. Crowding is also important during ri- 
bosomal translation on mRNA and binding protein 
diffusion along DNA 0, [1] , where bound proteins are hin- 
dered from passing each other. Furthermore, advances 
in nanofluidics allow studies of geometrically constrained 
nano-sized particles d, The system considered in 

this paper, the diffusion of a tagged particle immersed 
in a one-dimensional bath of hard-core interacting parti- 
cles - in the literature referred to as single-file diffusion 
(SFD) - represents one of the simplest systems governed 
by crowding effects, but yet with possible applications for 
obstructed one-dimensional protein diffusion along DNA 
molecules and transport in nanofluidic systems. 
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SFD phenomena emerge in quasi one-dimensional ge- 
ometries. The particle order is under these circumstances 
conserved over time t which results in interesting dynam- 
ics for a tagged particle, quite different from what is pre- 
dicted from classical diffusion (governed by Fick's law). 
Examples found in nature include ion or water trans- 
port through pores in biological membranes 

EH. 

one- 
dimensional hopping conductivity [12j | and channeling in 
zeolites [13| . SFD effects have also been studied in a num- 
ber of experimental setups such as colloidal systems and 
ring-like constructions [3, [H, [H, H3, 13 • One of the 
most apparent characteristics of SFD is that the mean 
square displacement (MSD) Sit) = ((yr — Ut,o) 2 ) 01 a 
tagged particle is in the long time limit proportional to 
i 1//2 in an infinite system with a fixed particle concentra- 
tion (brackets denote an average over initial positions and 
noise, yr and yr,o are tagged particle positions at time 
t and t = 0, respectively). Also, the conditional prob- 
ability density function for the tagged particle position 
(tPDF) is Gaussian. 

The first theoretical study showing the t 1 / 2 law of the 
MSD and that the tPDF is Gaussian is Ref. [H[. Subse- 
quent studies, proving the MSD law in alternative ways, 
are found in [2(| l2ll l22l [23l . Simple arguments to its 
origin are presented in [24l l25l |26| |. one of which [26| 
uses a simple relationship between the displacement of a 
single particle and particle density fluctuations, the lat- 
ter known to be the same as for independent particles 
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[l9l [27| . The i 1/,2 -law and Gaussian behavior have, in 
the long time limit, been shown to be of general validity 
for identical strongly over-damped particles interacting 
via any short range potential in which mutual passage is 
forbidden [28[. A generalized central limit theorem for 
tagged particle motion has also been proven [29(. More 
recent work include [30| where particles interacting via 
a screened Coulomb potential (i.e. not perfectly hard- 
core) was studied numerically (see also 131 1 1 . 133 d eals 
with SFD in an external potential, and [33l . |34| . 1351 ] ad- 
dress SFD dynamics with different diffusion constants. 
A phenomenological Langevin formulation of SFD was 
presented in pjrj ]. 

Although much work has been dedicated to single-file 
systems, to our knowledge, very few exact results for 
finite systems with finite-sized particles have been ob- 
tained. The one exception is [37| where the PDF for 
N point-particles diffusing on a finite one-dimensional 
line was derived. However, simplified expressions for the 
tPDF was only considered in the thermodynamic limit 
(A, L — > oo where L denotes system length and the con- 
centration g = N/L kept fixed). In this paper, we go 
beyond previous studies in the following ways: First, 
finite-sized particles are considered and we show that the 
iV-particle probability density function (APDF) can be 
written as a Bethe ansatz solution. We obtain an exact 
expression for the tPDF in terms of Jacobi polynomials 
(or hypergeometric functions) which reduces to that in 
|37| for the case of point particles. Second, we perform a 
(non-standard) large A-analysis of the tPDF, keeping the 
system size L finite. The expression for pr{VT, t\yr,o) m 
the many-particle limit is new and is presented compactly 
in terms of modified Bessel functions. An analysis of the 
tPDF reveals the existence of three dynamical regimes for 
a particle located roughly in the middle of the system: 

(A) short times, t -C r co ii = l/(g 2 D) where t co h denotes 
the collision time and D is the diffusion constant. In 
this limit the tagged particle undergoes standard Brow- 
nian motion with a MSD S(t) ~ t; (B) For intermediate 
times, T co ii <C t <C r cq where r cq = L 2 /D is the equi- 
librium time, we get a SFD regime where S(t) ~ t 1 / 2 . 
(C) for long times, t ^S> r oq , an equilibrium tPDF of 
polynomial-type is found. Notably, only regimes (A) and 

(B) exist in infinite systems. 

This paper has the following organization: Sec. HTlcon- 
tains the formulation of the problem and a mapping onto 
a point particle system. The tPDF is also formally stated 
in terms of the APDF to which governing dynamical 
equations are introduced. In Sec. IIIH we provide the 
solution to the equations of motion for the APDF using 
a coordinate Bethe ansatz. In Sec. IIVI the initial coor- 
dinates as well as the coordinates for all particles except 
the tagged one are integrated out in order to obtain an 
exact expression for the tPDF. Also, asymptotic results 
for large N for the tPDF are derived. In Sec. |V] the 
asymptotic large A expression for the tPDF is expanded 
for short and long times, and three different time regimes 
(A) — (C) (see above) are identified. More technical de- 
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FIG. 1: (color online) Cartoon of the problem considered here: 
N particles of linear size A diffusing is a one-dimensional sys- 
tem of length L. The particles have center-of-mass coordi- 
nates yj and initial positions yj.o, j = 1, ■•, N, and are unable 
to overtake. This implies that j/j+i > y-j + A (j = 1, .., N — 1) 
for all times. Also, the particles cannot diffuse out of the box, 
i.e. yi > -L/2 + A/2 and y N < L/2 - A/2. 

tails are given in the appendices. A brief summary of 
some of our results, corroborated with Gillespie simula- 
tions (Monte Carlo- type), can be found in [381 ]. 

II. PROBLEM DEFINITION 

In this paper we consider a system of A identical hard- 
core interacting particles, each with diffusion a constant 
D and a linear size A, diffusing in a finite one dimen- 
sional system extending from —L/2 to L/2. A schematic 
cartoon is depicted in Fig. [TJ The particles each have 
center-of-mass and initial coordinates y = (j/i, . . . , yjv) 
and yb = {yi,o, ■ ■ ■ ,Vn,o)i respectively. Due to the hard- 
core interaction, the particles cannot pass each other and 
retain their order at all times, i.e. Vj+i > yj + A, for 
j = 1 , . . . , N — 1 . The end of the system are reflecting 
(the particles cannot escape), i.e. y\ > —{L — A)/2 and 
Vn < (L — A)/2. 

The diffusion of finite-sized particles can be mapped 
onto a point-particle problem. Introducing the rescaled 
effective system length 

£ = L-NA, (1) 

and making the coordinate transformation 

•A 

X 3 = Vj-J A +— — A > 

Xjfl = %,o-jAH — A, (2) 

leads to 

K : - i/2 < Xi < x 2 . . . < x N < 1/2 (3) 

where 1Z denotes the phase space spanned by ([3]). The 
phase space TZq is also introduced for the initial coordi- 
nates which satisfy —i/2 < xi t o < £2,0 ■ ■ ■ < xn,o < £/2. 
For convenience, we also introduce the short-hand nota- 
tion x — (xi, . . . , xn) and x$ = (xi ; o, ■ ■ ■ , xn,o)- Equa- 
tions (QJ and ^ maps exactly the problem of TV finite- 
sized hard-core particles in a box of length L onto a N 
point-particle problem in a box of length £. 
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The main quantity of interest in this study is the tPDF 
Pt{xti t\xT,o)i that is the probability density that a 
tagged particle T (T = 1, N) is at position xt at time 
t, given that it was at xt,o at t = [an ensemble aver- 
age over over the initial (equilibrium) distribution of the 
surrounding N — 1 particles is implicit] . The equilibrium 
tPDF is straightforwardly calculated from the ergodicity 
principle: all points in the allowed phase space 1Z are 
equally probable. This leads the equilibrium iVPDF 



V eq (x) 



Nl 

JW 



N-l 



8(x i+ i - xt), 



where 9(z) is the Heaviside step function, 8(z) = 1 for 
z > and zero elsewhere. Using an extended phase space 
integration technique (see Appendix[C| it is easy to verify 
that f n dx\ ■ ■ ■ dxNV e(l {x) = 1. Integrating Eq. (|4]) over 
all coordinates leaving out one, xt, gives the tPDF 



dx\ ■ ■ ■ dx' N S(xr — x T )V eq (x') 



n 

Nl 1 



£ N NrlN, 



1/2 

x / dx'-r 



dx\ ■ ■■ dx' r _ x 

-112 J-i/2 



dx iv 



r+i / UUj n 

XT J XT 



Nr 



where S(z) is the Dirac delta-function and Nl {Nr) is 
the number of particles to the left (right) of the tagged 
particle (N — Nl + Nr + 1). In the remaining part of 
this section, we show how to calculate the complete time 
evolution of the tPDF from the many-particle iVPDF. 

In order to obtain pq-(xT, t\xT,o) one needs to first 
introduce the iV-particle joint probability density 
V{x, t; xq) which gives the probability density that the 
system is in a state x and that it initially was in a 
state xq. The joint probability density for the tagged 
particle pt(xt, t; a>r,o) is simply obtained from the joint 
iVPDF by integration over 1Z and TZq: pT{xT,t',XT,o) — 
J n dx[ ■ ■ ■ dx' N J no dx' lj0 • • • dx' N 6(xT - x' t )5(xt,o - 
x'q- )V(x', i; Xq). Relating the conditional and 
joint probability densities using Bayes' rule 
IH, i.e. V{x,t;x ) = V(x,t\x )V cq {x ) and 
PT(xT,t;x r ,o) = PT(xT,t\x r ,o)p cq ( x T,o), leads to 



pT{x T ,t\xT,o) 



1 



dxi ■ ■ ■ dx' N 6(xj- — x'q-) 



P oq 0*r,o) Jn 

dx[ o • • • dx' N S(xr : o - x' T ) 



xp(f,t|f )p^(4), 



(0) 



where p^(xr,o) is given in Eq. ([5]). 

In order to get the tPDF using Eq. ([6|), we need to 
calculate the iVPDF V(x,t\xo). It is governed by the 



diffusion equation 



dV(x,t\x ) 
dt 



D 



d 2 



d 2 



dx\ dx 2 



dx 2 N 



V{x,t\xa) 



(7) 

for x G 1Z (V(x, t\xo) = outside 1Z). The equation cer- 
tifying that neighboring particles cannot overtake reads 



D 



d 



i+l 



= 0. 



(8) 



(4) Also, reflecting boundaries are placed at the system ends 



D 



dV(x,t\x ) 



D 



dxi 



&P{x,t\x ) 



8xn 



x 1= -t/2 



x N =il2 



0. 



(9) 



(10) 



making sure that the particles are restricted to 
[— £/2,£/2] at all times. Finally, the initial condition is 



V{x, 0\x ) = 5(x 1 - x-Lfi) ■ ■ ■ S(x N - xn,o)- 



(11) 



Summarizing this section, the problem of ./V hard- 
core interacting particles of size A diffusing in a one- 
dimensional system of a finite length L was mapped onto 
a point-particle problem using relationships (fT]) and |2]). 
The dynamics of the ATTJF is governed by Eqs. ([7|-([I0|). 
Once they are solved (topic of Sec. IIII|) . the tPDF can 
be calculated via Eq. ([6]) (as demonstrated in Sec. IIV[) . 



III. iVPDF AS A COORDINATE BETHE 
ANSATZ 

In this section we obtain the iVPDF for the diffusion 
problem defined in previous section using a coordinate 
Bethe ansatz. The Bethe ansatz has been proven useful 
in solving a large variety of interacting particle problems 
since its introduction by Hans Bethe in 1931 (see [4(| 
for a review). The Bethe ansatz solution for the present 
problem reads 



V(x,t\x ) 



dkx [°° dk 2 



dk 



N 



2- 



x ^ e ik 1 x 1 +ik 2 X2+ik 3 X3 + ...+ik N x N 
_|_ g^ e ik 2 x 1 +ik 1 X2+ik; i X3 + ...+ik N x N 
I C Q c ik 2 x 1 +ik 3 x 2 +k 1 x 3 + ...+ik N x N 

+ all other permut. of {fci, k 2 , fc^v}] (12) 

where E(ki, . . . , fc^r) is the dispersion relation, 5V,- are the 
scattering coefficients, and 4>{kj,Xjfi) denotes a function 
containing boundary and initial conditions. Each one of 
these quantities are described below. 
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The dispersion relation has the form 



D(k\ 



+ k N ), 



(13) 



and relates "energy" to the momenta k\, . . . , fcjy- Equa- 
tion (fT3|) is obtained by inserting the Bethe ansatz Eq. 
([T2")) into the equation of motion, Eq. ([7]). 

The scattering coefficients Sy describe pair-wise par- 
ticle interactions and are in general functions of the mo- 
mentum variables fcj and kj, S!y = S(ki,kj). They are, 
however, independent of the initial positions of the par- 
ticles. In Appendix [A] it is demonstrated that the scat- 
tering coefficients making sure that the particles cannot 
pass each other [i.e. satisfying Eq. ((8])] are given by 



1, 



(14) 



which means that they are independent of momenta and 
correspond to perfect reflection. For non-interacting par- 
ticles Sij = 0, which reduces the Bethe ansatz to a stan- 
dard Fourier transform. 

The quantity (j)(kj,Xj t o) contains information about 
the initial and boundary conditions of the problem, here 
defined by Eqs. (HJ) - (fTTj) . The form of (p(kj,Xj } o) satis- 
fying these relationships is given by 



4>(kj , Xjfi) = 2 cos[fcj (xj t o + £/2)] 



oo 



^■(2171+1/2)1 



which is shown explicitly in Appendix[AJ Notably, for an 
infinite system we have <fi(kj,Xj. Q ) = e~ tkjXi -° |4ll. |43|. 

It is interesting to note that for SFD systems described 
by Eq. (fT2")) - (|15[) any macroscopic quantity that is in- 
variant under interchange of any two particle positions, 
Xi <-* Xj, takes the same value as for a system of non- 
interacting particles. This is in marked contrast to micro- 
scopic quantities such as the tPDF which in general be- 
have very differently for single-file and independent par- 
ticle systems. In Appendix IF1 we use the Bethe ansatz to 
explicitly calculate two macroscopic quantities, the dy- 
namic structure factor and the center of mass PDF, and 
show that they agree with standard results for indepen- 
dent particle systems. 

Integration over momenta in Eq. (fl~2|) [using Eqs. fl3])- 
([15])] leads to the iVPDF: 

V(x,t\x ) = ip(xi,xifl;t)ip(x2,x 2 fl;t)---ip(x N ,x N! o;t) 
+ ip(x 1 ,X2,o;t)ip(x 2 ,xi, Q ;t) ■ ■ ■ ip(x N , x N . ; t) 
+ all other permut. of {xi,o, %2,o, ■■■,Xn,o}, (16) 

where 

1 

ip(xi,Xj t0 ;t) 



(4 7 rDi) 1 /2 

OO s 

i exp 



exp 



(xi - Xjfi + 2ml)' 1 
(xi + x j}0 + (2m + l)i) 2 



4Dt 



,(17) 



is obtained from the inverse Fourier transform 
(27r)- 1 /^° oo dfc J -^{k^Xj.^e-^e 1 ^. We point out 
that ip(xi,Xj } o; t) is the single-particle PDF for a particle 
in confined in a box of length I. 

The single-particle PDF given in Eq. (jTTJ) is, however, 
not convenient for analyzing the long time limit t — > oo. 
In order to get a more suitable expression we seek instead 
the eigenmode expansion of ip(xi, Xj^;t), which can be 
done in a variety of ways. Here we use Bromwich in- 
tegration. The Laplace transform of Eq. (fT7|) is (see 
Appendix iBj) 



ip(xi,Xj }0 ;s) 



die st ip(xi,Xjfl;t) 



V4Dssinh(e^/s/D) 
x ^cosh[(xi + Xj,o)\/ s/D] 

+ cosh[(£ - \x jfi - x^^fl/D]) .(18) 

The sought eigenvalue expansion is obtained as a sum of 
residues of ip(xi, Xj : o; s) (see e.g. Ref. HU), and reads 



ip(xi,Xj j0 ;i) 
where 



1 



1+5^ Gm(Xi,Xjfi)£m{t) (19) 



m— 1 



r ( \ (+) frnnx t \ / 

y m (xi 7 x J} o) = v^'cosy — - — j cos ^ 



mTTXjfl 

£ 

mirxjfi 



i \ . /mirxi\ . /rmTXjo\ ,„„ N 



2 Dt/i 2 



and 



,(±) 



1±(-1)' 



(21) 



(22) 



Elementary trigonometric identities [44J were used to 
bring Q m {xi, Xj.o) onto the form in Eq. (I20p . Equations 
(fT9| - ([22|) agrees with well-known results [4JJ. The single- 
particle PDF (jTSJ) is more convenient for obtaining the 
long time limit as well as for numerical computations 
compared to Eq. (fT7|) . 

In summary, the many-particle APDF for excluding 
particles of size A diffusing in a finite interval of length 
L with reflecting boundaries is given by Eqs. (TT^j) . HH) 
[or Eq. (jTTJ)] combined with the mapping equations (JTJ) 
and (|2|). For point-particles (A = 0) these results agree 
with those presented in [37| where a different approach 
was used [46j ]. Based on the explicit expression of our 
A^PDF, we will in the following section address the tPDF. 



IV. tPDF - EXACT AND LARGE N RESULTS 

In this section, we calculate the tPDF © by inte- 
grating out the coordinates and initial positions of all 
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non-tagged particles from the 7VPDF given in Eq. (|T6|) 
[47l ]. As is shown in detail in Appendix [Q we can, due 
to the property that V(x,t\xo) is invariant under per- 
mutations of X% 4 — ^ Xj -j extend the integration from 1Z 
to the hypercubes Xj £ [—£/2,xt], j — 1 and 

xj 6 [xt,£/2], j = T + 1,...,N. A similar procedure 
holds for integration over TZq (initial positions). Using 
the extended phase-space technique, Eq. ^ becomes 



PT(xT,t\x T fi) 
r £/2 



f N L f N R 
J L JR. 

NrlNpl 



R- J- 



dx\ 



d,XT- 



1/2 



dx 



N 



■I'T 



XT,0 



dxifi ■ 



(23) 



\ I dxr+i 

' XT 

\ / dxT-1,0 dxr+1,0--- dx Nfi V(x,t\x ), 

-1/2 Jxt,o «I,0 



where 

f L = (e/2 + x T , )-\ 



f R = (1/2 - x T ,oY 



(24) 



Using a similar combinatorial analysis to the one in Ref. 
[53], we arrive at Eq. (iDlj) (see Appendix H3J. The 
tPDF is, however, more conveniently expressed in terms 
of Jacobi polynomials (44[, P^^\z). Using the identi- 



ties P, 



(z) = (n + a)\(n + /3)\/[n\(n + a + 0)\}[(z 



l)/2]- a P}ig p >(z) and P^ p) (-z) = (-l)"P^ a >( z ) [4 
leads to 

p(xr,t\xr,o) = 

x{(N 



iNL + NR - 1)l ^ L ) N ^ R R.) N « 



N L IN R \ 



n 



$(1,0, 0;0 



+N; 



■*(0,l,0;O 



-N l Nr 



(0,0, 



(25) 



where 



$(a,6,c;0 = 



(N L -(a + c)y.{N R -b)\ 



(N L + N R -(a + b+c))\ 

X £ C (l - C)N L -(a+c) p(c,N R -N L +a-b) (\ 



-i 



and 



(26) 



(27) 



were introduced. The quantities tpL, ip^ etc. are defined 
in Eq. (|D3|) and are integrals of the one particle propa- 
gator ip = ipi xj, Xjfi, t) [43]. For the general case we have 
i <E [0,1] [50]. By using limiting results of ip^, ip R , ip^ 
and tpR from Appendix [El one concludes that £ — > for 
short times it — > 0), and £ — » 1 in the long time limit 
(t — > 00). It is also possible using standard relations for 



the Jacobi polynomials [44| to express <E>(a, 6,c;£) as a 
Gauss hypergeometric function 2 Fi {ol, P, 7! z): 

$(a,6,c;£) = 2 Fi(-N L + a,-N R + b, 

~(N L +N R )+a + b + c;l-0 (28) 

which is convenient for obtaining the long-time behavior 
as will be seen in the next section. 

Finally, we need to find explicit expressions for the 
integrals of ip(xi, Xj t o] t) [Eq. (|E)3|) ] . Integrating Eq. ([T9[) 
[or integrating ip(xi, Xj^; s) prior to Laplace inversion, 
see Appendix [E] yields (arguments left implicit) 



1 XT 

2 + T 



1 xr.o 

2 1 



i>R 

i> 

IpL = 

IpR = 
V* = 

with 



1 XT (1 XT.Q 

2 £ V 2 £ 



m— 1 

— 1 00 



1 °° 



1 _ ^r,o 
2 



2J J m (xT,Q,XT) £m(t) 



I 



111 — 1 

— 1 OO 



J m {xr^xr)S m {t) 



711 — 1 



1-^,^ = 1-^,^ = 1"^ (29) 



1 



(run) 2 



-vL ' cos 



f+ 1 . (mTTXT\ . (m-KXTfi 
/mirxT\ I 



( ^ ) ^ ) 



1 

TO7T 



i^i^ t ~- ) sin 



/muz 



— i>L ' cos 



/mirz\ 



rmtXTfi 
/ mnz' 



cos 



TO7TZ 



(30) 



(31) 



where £ m (t) and ^^f^ are given by Eqs. ([2"Tj) and 
respectively. To summarize, the complete expression for 
the tagged particle PDF is given by Eqs. p^[) -(|27 [l and 
([2^ ]) -(|5D |l together with Eqs. © and © for the case of 
finite-sized particles. The expressions for pt(xt ,t\xT,o) 
can straightforwardly be computed numerically; our 
MATLAB implementation is available upon request. 

In the remaining part of this section we derive the 
tagged particle PDF Pt{xt ,t\xT,o) valid for a large iV" 
and (finite) system size I. This large N expansion will 
be used in the next section for identifying different time 
regimes and to obtain pt(xt, t\xT,o) f° r short and inter- 
mediate times. From Eq. ([2"6"[) we note that the argu- 
ment in the Jacobi polynomial, PlC'^iz), is in the inter- 
val z £ [l,oo) (since £ £ [0, 1]) and that the number of 
particles is related to the order n. A large TV expansion 
of $(a, 6, c;£) therefore amounts to find a large order n 
expansion valid for z £ [l,oo) (i.e. for all times) for the 
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Jacobi polynomial. One such expansion was derived in 
[5l| (see also [H|,[53j]), and applying it to Eq. (j2"6"j) yields 

$M,c;0 « [N -{a + b + c)] 1 ' 2 ^ 2 ^^ 



x(2ttC) 



1/2 



[JV-(a+6+c)]/2 



x/ c [(iV-(a + 6 + C ))C](l + A(C)) 



where [54 



c 



111 



i + Vf 



(32) 



(33) 



and J Q (.z) is the modified Bessel function of the first kind 
of order a. Stirling's formula [44| was also used to ap- 
proximate factorials involving Nl and Nr. The correc- 
tion term appearing in Eq. (|32p is 



MO 

Bo(0 



Bo(0 



l c+l 



[(N-(a + b + c))(] 



N-(a + b + c) I C [(N- (a + b + c))(] 
i{(c 2 -l/4) (j-cothc) 

-[(fo + a-&) 3 -l/4]tanh(C)}, (34) 



where Sn = Nr — Nl was introduced. The expression 
for A(Q was found by explicitly evaluating an integral in 
Ref. [5l| , and it is a straightforward matter to show that 
the correction term A(() is indeed always small for all 
£ € [0,1] provided that 1/N, \S N \/N < 1. Inserting Eq. 
(|32jl in Eq. (|25p and using the same approximations as 
above, we obtain our final large N result for the tPDF: 



W>l ) ( n ' Sn ~ 1)/2 (iPr) {n+5n ~ 1)/2 

1/2 



PT(xT,t\xT,o) 

x(i-o (JV - 1)/2 (^) 1/z {(i-e) v Vo[ivc] 



+ 



iV 



lj} R 



v 7 ! 



Io[N( 

R 

V^R 

i> L R 



h[NQ). 



(35) 



We point out that this expression, in contrast to pre- 
vious asymptotic expressions [2(J H3] , is valid for a 
finite box of size £, assuming only that the number of 
particles is large f > 1 and that the tagged particle is 
approximately in the center of the system: \5n\/N -C 1. 



V. THREE DIFFERENT TIME REGIMES 

In this section we show that for large N, the finite SFD 
system considered here has three different time regimes to 
which expressions for pq-(xr, t|a?T,o) are derived. Math- 
ematically, the different cases appear due to the magni- 
tude of N( [found in the argument of the Bessel functions 



in Eq. (|35[)] . and if ( is small or large. Utilizing Eqs. 
(|E14[) and (j3"3"|) . these cases can be turned into different 
time regimes if introducing the collision time 



Tcoll 



Q 2 D' 



(36) 



where g = N/£ is the concentration of particles, and the 
equilibrium time 



D ' 



(37) 



For a particle located roughly in the middle, \8n\/N <C 1 
(i.e. Eq. (|35j) applies), the three cases are given by 

A. short times, N( <C 1, i.e. t <C t co h, r oq , 

B. intermediate times, £ <C 1 and iV£ 3> 1, corre- 
sponding tO T co n <C t < T C q, 

C. long times, C 3> 1, i-e. t r co n,T cq . 

Time regimes (A)-(C) are analyzed in detail below. 

A. Short times, t <§C Tcoli, ^~cq 

For short times we have N(^ <C 1, and may therefore 
use the approximations / a (z)| z <i ~ (z/2) Q /r(a+l) [H ]. 
C ~ VS, an d ^ < f In this limit, one finds that the first 
term in Eq. (|35[) dominates which in combination with 
Eq. © leads to 



Pt(ut jt\yT,o) — (AirDt) 1/2 exp 
for which the MSD is 

S(t) = 2Dt. 



(yr ~ yr,of 
ADt 



(38) 



(39) 



In the short time regime, almost no collisions with the 
neighboring particles (nor the box walls) have occurred 
and the tPDF is therefore a Gaussian with width 2Dt as 
for a free particle in an infinite one-dimensional system. 



B. 



Intermediate times, r co ii < t < t ci 



In the intermediate time regime the tagged particle 
has collided many times with its neighbors but not yet 
reached its equilibrium tPDF. For this regime, where 
C < 1 but N( > 1, we get the tPDF as follows. First, 
the first term in Eq. (|3"5"|) is neglected (this is checked 
at the end of the calculation). Second, the Bessel func- 
tion is approximated with I a {z)\ z ^i ~ e z /y/2nz. A 
straightforward expansion of Eq. (|3"5"|) for ip^, ip^ <§; 1 
(i.e. v 7 ? 1), together with Stirling's formula, gives: 



PT 



7 



^ WL + l/> R ll>R 



(40) 



If we furthermore assume t hat t he average of the absolute 
value of 7] — (xt —xt,o)/ %/ &Dt is small (which is checked 
after the calculation), Eq. (|E15[) may be used which in 
combination with Eqs. pO]) , |T]) and ([2]), keeping only 
lowest order terms in 77, leads to the SFD result: 



1/4 



PT(yT,t\yr,o) 



x exp 



where the MSD is 



(VT - VTS)) 2 



(41) 



1 - gA 4Dt 



S(t) = — ^\ — 



(42) 



Equation ([42|l justifies the assumption that the expec- 
tation value of |?7| is a small number. Also, comparing 
the magnitude of the first term in Eq. (|3"5")1 with respect 
to the second and the third, shows indeed that our first 
assumption above was correct 55|. For point-particles 
A = 0, Eq. (|4"Tj) agrees with standard results [20, H3|) m 
which N, L — > 00 while keeping the concentration g fixed. 
The result above shows that SFD behavior appears also 
in a finite system with reflecting ends, as an interme- 
diate regime (for a particle roughly in the middle). In 
addition, note that the simple rescaling g — » g/(l — gA) 
takes us from previous point-particle results to those of 
finite-sized particles. 



C. Long times, t 3> T cq 

In the long-time limit we have £ 3> 1, for which the 
tPDF can be obtained exactly for arbitrary N. Using 
the exact expression for p<r{xr, t\xr,a) found in Eqs. (|25)) 
and pgjl. together with 2 F 1 (a, (3,j; z = 0) = 1 [4^ and 
Eq. jEJJl [also using Eqs. (JTJ and ©] gives 



Pt(vt) 



1 



(n l + n r + i)\ 



(L-NA) N NlINrI 



x [ - + y T - A(l/2 + N L ) 



L 



x ( - - y T - A(l/2 + N R ) ) . (43) 



where p^{yr) = PT(yT,t — * oo\yT,o)- This equation 
agrees with the equilibrium tPDF given in Eq. (J5j) (for 
A = 0), as it should [56[. The results above shows the 
consistency of our analysis, and illustrates the ergodicity 



of the finite SFD system. Calculating the second moment 
of the equilibrium tPDF, S(t — * 00) = <S oq , gives 



S e q — 



L- NAY 



r(i/2)r(2(jv fl + 1)) 



T{N R + 1)T(N R 



5/2)' 
(44) 

for Nl = N R where T(z) is the the gamma function. For 
N 3> 1 we can simplify the expression for the equilibrium 
tPDF as well as S cq . If a symmetric file is assumed, i.e. 
Ni = N R and T = (N + l)/2, an asymptotic expansion 
of p^(yr) (using Stirling's approximation [44|) gives the 
Gaussian PDF 



<2N 



1 



7T (L — NA) 2 



x exp -2N 



L — NA 



from which we read off that 

(L - NA) 2 



AN 



(45) 



(46) 



Equation (|46|) can also be found directly from a large 
N expansion of Eq. (|44[) using Stirling's formula and 
assuming N 27Vr. We point out that p^(yr) — ► for 
N, L — > 00 even if the concentration g = N/ L is kept 
fixed. This is consistent with the long time limit of Eq. 
(|4"T|l which indeed goes to zero for large times. Finally, 
the analysis in this subsection also gives an estimate on 
the time required to reach equilibrium, namely t 3> T oq . 



VI. CONCLUSIONS AND OUTLOOK 

In this study we have solved exactly a non-equilibrium 
statistical mechanics problem: diffusion of N hard-core 
interacting particles of size A which are unable to pass 
each other in a one-dimensional system of length L 
with reflecting boundaries. In particular, we obtained 
an exact expression for the probability density function 
Pr(yT,t\yT,o) (denoted tPDF) that a tagged particle 
particle T is at position yj- at time t given that it at 
time t = was at position yr,o- We derived the tPDF by 
first finding the iV-particle probability density function 
(iVPDF) via the Bethe ansatz, and then integrating out 
the coordinates and taking the average over the initial 
positions of all particles except one. The exact expres- 
sion for pr{yT,t\yTfi) is found in Eqs, fT}, (f2"S"]) and 
(|26|) . and constitutes the main result of the paper. For 
a large number of particles and for a tagged particle lo- 
cated roughly in the middle of the system, an asymptotic 
expansion of the tPDF was derived [see Eq. ([55])]. Based 
on this equation, we found three time regimes of interest: 
(A) For short times, i.e. times much smaller than the col- 
lision time t <C Tcoii = 1/(q 2 D), where g = N/L is the 
particle number concentration, the tPDF coincides with 
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the Gaussian probability density function which charac- 
terizes a free particle [Eq. (B) For intermedi- 
ate times, t 3> T co n but much smaller than the equilib- 
rium time t <C r cq = L 2 /D, a sub-diffusive single- file 
regime was found in which the tPDF is a Gaussian with 
an associated MSD proportional to t 1 ! 2 [Eq. (C) 
For times exceeding the equilibrium time t 3> T eq , the 
tPDF approaches a probability density of polynomial- 
type [Eq. (J43J)]. 

We point out that the sub-diffusive behavior for a 
tagged particle in time regime (B) is of fractional brow- 
nian motion type [33|, H3, rather than that of 
continuous-time random walks (CTRWs) characterized 
by heavy-tailed waiting time densities [53 . loll l6l| . For 
such CTRW processes the probability density function is 
not a Gaussian as for the current system. Further com- 
parisons between sub-diffusion in single-file systems and 
that occurring in CTRW theory was pursued numerically 
in [El. 

The Bethe ansatz is often employed in quantum me- 
chanics when many- body systems are studied (e.g. quan- 



APPENDIX A: THE BETHE ANSATZ 

In this section we show that the Bethe ansatz, Eq. 
(fT2")l . is a solution to the problem defined by Eqs. ([7])- 
(fTTj) . First, it is demonstrated that Eq. (fT2"| satisfies the 
boundary conditions at the ends of the box. Second, we 
show that the requirement that the particles are unable 
to pass each other is satisfied by setting the scattering 
coefficients to unity. Finally, it is demonstrated that Eq. 
(fT2l) also satisfies the initial condition. 



1. Boundary conditions at the ends of the box 

In this subsection it is proven that Eq. (fT2"|) satisfies 
the the reflecting boundary conditions Eq. ^ and (fTU|) 
at ±£/2 with an appropriate choice for <p(kj,Xjfi) [Eq. 
(fl~5|) ] . The scattering coefficients are set to SV,- = I which 
is proven to be correct in the following subsection. First, 
we define the function 



turn spin chains) [40|, l63Lanch also for stochastic many- X (k, z) = 6(k, z)e" iM/2 = 2co8\k(z+£/2)] V 

„„„-t;„i„ i„j-j-;„„ — „ui„„,„ Ic/i cell \\r„ u r j-T,„± j-i,„ +^„„ ^ — ' 



particle lattice problems [6J, [65( . We hope that the theo- 
retical analysis based on the Bethe ansatz presented here 
will stimulate further progress in the field of single-file 
diffusion and that of interacting random walkers. For 
instance, it would be interesting to see whether our anal- 
ysis could be extended to derive exact results also for 
particles interacting through potentials other than hard- 
core type l3l| . and for particles having different diffusion 
constants (34j . 

From the applied point of view, our exact expres- 
sion for the tPDF covers all time regimes and is 
straightforward to implement for numerical computa- 
tions. Therefore, we believe that our explicit formula 
for pt(uTi t\VT,o)i as well as the approximate results for 
regimes (A) — (C) > wn l be useful for experimentalists (see 
for instance [18|) seeking to extract system parameters 
such as particle size A, system size L, the particle's dif- 
fusion constants (D), and the number (Nl and Nr) of 
particles to the left and right of the tagged particle. 

We finally note that the use of fluorescently labeled 
(tagged) particles is of much use for studying biological 
systems. The understanding how the motion of such par- 
ticles correlate with its environment is therefore the key 
for grasping the behavior of such systems in a quantita- 
tive fashion. 
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which has the symmetry relation 

X(k, z) = A(— k, z). 



~2ikrnl 



(Al) 



(A2) 



By taking the derivative of Eq. (fT2")) with respect to ii, 
and evaluating at the left boundary x\ = —£/2 gives 



dP(x,t\x ) 



dx\ 
3 dk\ 



x 1= -l/2 

dk 



2ir 



ik ie iklXl (e ik '~ 



X2 + ... + ik.NX N 



^ = l4>(kj,Xj fi ) 

perm, of {fc 2 , k N }) 



+ ik e ik ' Xl ^ e ik 2X2 + ---+ik 1 x i + ...+ik N x N 

+ perm, of {k 2 , fej-i, k i+1 , k N }) 



ik^e ikNXl (g ik ' i X2+---+ikix I , 



perm, of few-i}) 



x 1= -t/2 



dk 2 
"2V 



dk 



2n 



f ^(iki)e- D ^X(k 1 ,x lt0 ) + 



+ 



dk\ 
~2tt 



dki-i f°° dk 



-00 

e -u(fc?+...+fc?_ l +*? + x+...+^)*njL 1 , ¥i 0(fc 



2lT 



2n 



00 dk jv 



2tt 



3 = l^^iV\^3 1 -^3 
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x r P(*k l )e- Dk ? t \(k l7 x l ,„) + 

J-oo 27T 

J-oo 2tt 7„ 00 2ir 
xn ] t 1 Vft 1 i j ,o) ^ ^-(ik N )e- Dk ^X(k N ,x N , ) 
= 0, (A3) 

where f^° oci dkikie~ Dkit \(ki,Xi.o) — [odd integrand, 
see Eq. (IA2[) ] was used in the last step. Note that 
this calculation does not rely on any specific form of 
(j)(kj,Xjfi). It is only required that the symmetry re- 
lation (|A2|) holds. Furthermore, the dispersion relation 
E(k) = D(k\ + . . . + k 2 N ) was also used in the above 
derivation. However, it would work equally well for any 



dispersion relation as long as E(k) = J2i E(ki) with 
E{ki) = E{-h) is valid. 

A similar analysis as just presented shows that the 
Bethe ansatz solution also satisfies the reflecting condi- 
tion at +1/2 [Eq. (fTU)) ]. In fact, since the Bethe ansatz is 
invariant under the coordinate transformation 
gives [dP{x, t\xo)/dXj] x .—±i/2 = for all Xj. 



2. Single-file condition: particles are unable to 
overtake 

In this subsection it is shown that the condition that 
the particles are unable to pass each other, Eq. ([5]), is 
satisfied for scattering coefficients given by SV, = 1 in 
the Bethe ansatz solution Eq. (fT2")l . We start off by 
expressing V(x, t\xo) in two alternative ways: 



V(x,t\x ) = 





00 dki 


I 


oo 27r 


X 


iki'Xj 


+ 





dki 
~2^ 



dk 



N 



2tt 



e- D «+- +fc -)*nf = ^(fe i ,^ ) 



all other permut. of {k 2 , /cat}) 



and 

V{x,t\x ) 



+ e ik N x ^k 1 x 1 + ...+ l k^ 1 x^ 1+l k ] + 1 x ] + 1 +ik J+2 x ]+2 + ...+ik ] x N + aU othcr permut of ^ k N -l}) (A4) 

r°° dh r°° dk2 r°° dkx D{k *+... +k %)t U N m . ) 

J-oo 2vr ./ x 2vr J x 2n 6 H^xm^o) 

e ik lXj+1 ^k J + 1 x 1 + ...+ t k^ 1 x^ 1 + t k J x ]+i k J+2 x ] + 2 + ...+ l k N x N + aU other permut of {k 2 ,...,k N }) 

1 x^ 1 +ik J x ]+l k ] + 2 x J+2 + ...+tk N x N + aU Qther permut of fc;Jj ; fcjv ^ 



_|_ e ik 2 x j+1 ^ 



^ik\x\-\-ikjj t \x 2 -\-...-\-ikj — 



+ ^k N x ] + 1 ^ l k 1 x 1 + ...+ l k^ 1 x^ 1 + l k ] x ] +rk ]+2 x ] + 2 + ...+ l k 1 + 1 x N +aU othcr permut of {k u ...,k N -!}) . (A5) 

Using the above equations one finds 



/ dV(x,t\x ) dV(x,t\x ) 



V dxj+i 



dx; 



Xj + i —Xj 



dki 
~2^ 



dk 



E e -D(kl...,k%)t n N (f)ikjiXjo) 



-oo 2lT J-oo 27r 
ikie iklXj (^ l o i k ] + ixi + ...+ikj-iXj-i+ikjXj + i+ikj+ 2 Xj+ 2 + ...+ikNXN 
_ e ik ] x 1 +...+ l k ] _ 1 x ] _ 1 + l k ] + 1 x ] + 1 + l k ]+2 x ]+2 + ...+ l k N x N + aU Qther permut of ^ f, N ^ 

+ ik2e ik2Xj ^g iklXl + ik i+ lX2 +---+ ik 3- lX 3- 1 + ik 3 x 3+ 1 + ik 3+ 2X j+ 2 +---+ ikNXN 

_ e ik 1 x 1 +ik J x 2 + ...+ i k^ 1 x^ 1+i k J + 1 x J + 1 +ik J + 2 x J+2 + ...+ l k N x N + aU Qthcr permut of fc 3 , k N }) 



_ e ik 1 x 1 + ...+tk J - 1 x J - 1 +tk J + 1 x J + 1+i k ]+2 x J + 2 + ...+ik J x N + aU other permut of { fcl; ..^ jfc JV _ 1 j.) ] = 

I 



(A6) 



where it was used that each parenthesis is identically zero tation over all allowed momenta. Note that this deriva- 
due to the cancellation of the 2(^—1)! terms after permu- 
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tion is independent of the choice of <j)(kj, x^o) and E{k). 



3. Initial condition 

In this subsection we show that the Bethe ansatz Eq. 
(fT2|) agrees with the initial condition Eq. (jTTJ) when 
t — ► 0. By defining 



(A7) 



Eq. (H2J) reads 

V(x,t — ► 0|x ) = Xi )0 )fi(x2, 2:2,0) ■■ • ^(aJjv, Xjv )0 ) 

+ n(xi, 2:2,0)^(2:2,^1,0) • • • si(x N ,x Nt0 ) 

+ all other permut. of {xifi, xjv.o}-(A8) 



Using the explicit expression for 4>(k, y) found in Eq. |T5 
and that 5{x - z) = (27T)- 1 J™ dk e^C*-*) leads to 



fl(x,y) = 5(x-y + 2m£)+S(x + y+(2m+l)£). 

7TI — — OO 

(A9) 

For all to 7^ the ^-functions are non-zero for coordi- 
nates lying outside of 1Z and TZq where, per definition, 
V(x,t\xo) = (see discussion in section Sect. HI]) . For 
to = 0, it is only the first 8- function in the sum, 5(x—y), 
that contributes to the iVPDF. Furthermore, since all 
terms except the first one in Eq. (|A8[) are zero due to 
the fact that V(x,t\xo) — outside 1Z, one obtains 



V(x, t -* 0|x ) = S(x - x ) • • ■ S(x N - x Nfi ) 
as t — > which is the desired result. 



(A10) 



APPENDIX B: LAPLACE TRANSFORM OF 

i>(xi, x^o; t) 

In this section it is shown explicitly how one can go 
from the one-particle PDF ip(xi, Xjfi', t) for a particle in 
a box expressed in terms of Gaussians [Eq. (jTTJ)] , to the 
eigenmode expansion [Eqs. (T!?)) - (f2"2"]) ] used in this paper. 
First, the summation in Eq. (|17[) is divided schemat- 
ically into Em=-oo/™ = /m=Q + Em=l[/m + f-m]- 
Then, using the Laplace transform £[(7rf) _1 / 2 e~ a /( 4 *)] = 
s -i/2 e - Qv /i ( a > 0) [|3 one finds 



ip(xi,Xj fi ;s) = Q(s) + 



"(Si— Xj,o)v/ s/D , (xi-Xj^y/s/Dl 



,-(x t +x jl0 +&)y/<>/D _|_ e {x t +x j:0 +e)^/s/D 



x J2 e 

m— 1 



-ImlJ sjD 



(Bl) 



X, _ x 2 




FIG. 2: Integration area (the darker upper area, x\ < X2) 
for the integral defined in Eq. (|C2|I . For any function which 
is symmetric under «-» £2 the integration over the lower 
triangle (xi > £2) yields the same result as integration over 
the same function over the upper triangle (xi < X2). One 
may therefore extend the integration area to the full rectangle 
above provided one divides the corresponding result by 2. 



where 

Q(s) 



I (4Ds)- 1 /2 e -(^ l -^,o)V s /- D , Xl > x jfi 
\ (4L>s)- 1 /2 e -(^,o-^)V I 7 : D, Xl < x jfi . 

Considering the cases Xi > Xj t o and Xi < Xj t o separately, 
and that £„=i e - 2mt ^ = ( e 2 V^- l)-i , leads to 



(B2) 



sDsmh[£y/s/D] 

X COSh[(^/2 — Xj$)y/ s/D\, Xi < XjQ 



x < 



coah.[(£/2-Xi)y/s/B\ 
x cosh[(£/2 + xjfi)y/s/D], Xi > Xjfl 



(B3) 



which after elementary trigonometric manipulations re- 
sult in Eq. flgp . 



APPENDIX C: EXTENDED PHASE SPACE 
INTEGRATION 

When the tPDF is integrated out from the iVPDF we 
need to resolve the following type of integral: 



n\x-i 



T(xt) = \ dx[ ■ ■ ■ d X ' T _ 1 dx' T+1 ■ ■ ■ dx' N V(x' ,t\xo) 

(CI) 

over the region TZ [Eq. ©] with the tagged particle co- 
ordinate xt left out. As pointed previously in the pa- 
per, the Bethe ansatz solution [Eq. Q12p] is symmetric 
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under the transformation when all scattering 

coefficients are given by = 1. This allows the in- 
tegration of V(x,t\xo) in Eq. (|C1|) to be extended to 
the whole hyperspace Xj G \—l/2,xr\,j = 1,...,T — 1 
and Xj £ [xT,i/2],j = T + 1,...,N. This is most easily 
demonstrated in an example which then is extended to 
the general situation. Consider the case of three particles 
where particle three is tagged T = 3: 

Z{xg,) = / dxidx 2 V(xi,X2,x 3 ,t\x ). 

J -£/2<x 1 <x 2 <x 3 <£/2 

The integration area in the (x±, a;2)-plane is sketched in 
Fig. [2] (upper dark triangle). Since V(x,t\xo) is invariant 
under Xi «-> X2, integration over the lower triangle gives 
the same result, i.e.: 

Z( x a) = / dxidx2V{xi,X2,x 3 ,t\x ). 

J -£/2<x 2 <x 1 <x 3 <l/2 

(C3) 

If Eqs. (|C2[) and (|C3|) are added, T(x 3 ) can be expressed 
as an integral over the full rectangle 

X(cc 3 ) = — / dx 1 / dx 2 V(x l7 X2,x 3 ,t\x ) (C4) 
z J-e/2 J —1/2 

For the general case, the integration can be extended for 
any particle number to the left of the tagged particle. 
This means that it is possible to go from integration over 
the phase space —lj2<x\ < ... < xt-x < xr to — £/2 < 
Xj < xt for j = 1, . .. , T — 1 provided that we divide by 
Nl\. This holds also for Nr particles to the right and 
Eq. (jCip can in general be written as 



l(x T ) 



1 



xt r x T 

dx\ j dx 2 ■ ■ ■ (C5) 



N L \N R \ J_ i/2 J_ e/2 
fir rl/2 I- 1 / 2 

x / dxr-i / dxT+f" / dx^V{x, t\x*o). 

J-e/2 Jx T Jxt 

Since V(x, t\xo) is also invariant under Xi : o <-> Xj.o, a 
similar extended phase-space technique is valid for inte- 
grations over the initial particle positions TZq. 



APPENDIX D: FROM THE iVPDF TO THE tPDF 

In this section the tPDF is obtained from the integral 
([6]) using the Bethe ansatz iVPDF (p~6|) explicitly. A sim- 
ilar combinatorial analysis to that found in |37[ gives 



min{N L ,N R -l} 
9=0 

min{N L -l,N R -l} 

+ Yl H*{q){rt) N --«-\i, L R y+^ L 

9=0 

min{N L -l,N R -l} 
9=0 



(Dl) 



with the combinatorial factors 



N L \N R \N L 



and integrals (leaving arguments xt,o and xt implicit) 

/XT r x T,o 
dxi / dxjfli>(xi,Xj >0 ;t) 
-i/2 J-e/2 
r e/2 .1/2 
= fn / dxi / dxj t0 ip(xi,Xj t0 ;t) 

Jxt «t,o 
rZ/2 rxT.a 

V'lW = II dx % / dxj fi ip(xi,x jt0 ;t) 

Jxt J-e/2 
rx T j-e/2 

^ItCO = /h / / dxj t ot/j(xi,Xjfi;t) 

J-e/2 Jx T ,o 

/XT 
dxitp(xi 7 Xj fi ;t) 
-i/2 

ni/2 

ip R (t) = / dXitp(xi,Xj i0 ,t) 

J XT 

/XT,0 
dx j>0 ip(xj,Xj fi ;t) 
-i/2 

r e/2 

V'flW = /.R / dxjfi^(xj,x jfi ;t) 

Jxt,0 



(D3) 



minliVi.iVB} 

]T ^(g)^)^" 9 ^) 9 ^^) 9 ^)^- 9 

9=0 

min{iV£ — l,iVii} 
9=0 



The prefactors /r, and are found in Eq. and 
correspond to uniform distributions to the left and right 
of the tagged particle according to which the surrounding 
particles are initially placed. They appear when integrals 
over initial coordinates are performed. Also, it is easy to 
see from normalization that i/j L (t) + i/j R (t) = 1, ip^(t) + 
$ft) = land$[(i)+^(i) = :L: 
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If considering a single particle in a box of length £, 
the integrals defined in Eq. (|D3[) are easily interpreted 
as follows. First, ipj^ (i/jr) is the probability that single 
particle is to the left (right) of xt at time t given that it 
started, with an equal probability, anywhere to the left 
(right) of xt.o- Similar interpretations hold also for tp^ 
and i/jr. The quantity i/jl (V'i?) is the probability that 
a single particle is at position xt given that the particle 
started somewhere to the left (right) of xt,o- Finally, 
(^/,-R) j s probability that a single particle is to 
the left (right) of xt at time t given that it started at 
position xt o- 



APPENDIX E: INTEGRALS OF ^(x^x^a) IN 
LAPLACE-SPACE 



In this section, exact expressions as well as limiting 
forms for the integrals appearing in Eq. (|D3j) are given 
in the Laplace-domain. Using the Laplace-transformed 
one-particle PDF ip(%i,Xj,o',s) found in Eq. (|F33|) . the 
integrals defined in Eq. (|D3|) (with time t replaced by 
Laplace variable s, and x-r and xt,o left out) are given by 



1 



ssiDh[£y/a/D] 



' smh[(e/2 + x T )y r s/D} 

x cosh[(£/2 — XT,o)ys/D], xt < ir,o 



x < 



sinh[f y/s/D] - sinh[(£/2 - x T ) y/s/D] 
x cosh[(£/2 + xr,o)y/ s/D], xt > it,o 



(El) 



h 



!sinh[iV 's/D] 



' sinh[£ y/s/D] - cosh[(£/2 + xt) y/s/D] 
x sinh[(£/2 — XT,a)y/s/D], xt<xt.q 



x < 



^r(s) 



osh[(£/2-x T ) y/s/D] 
x sinh[(£/2 + xt,o)\/s/D], x T >x T .o 



frt 



(E2) 



s sm 



h[£yfi/D] 



' cosh[{£/2 + x T ) y/s/D] 



x sinh[(^/2 - x T ,o) y/s/D] , x T < x T ,o 



x < 



smh[£*/s/D] - cosh[(£/2 - xt) y/s/D] 
t x sinh[(£/2 + x T ,a) y/s/D], x T > x T ,o 



(E3) 



ri(s) 



s smh[£ y/s/D] 



f (1 - (xt,o - XT)h) smh[e x /7/D] 
-f L J%ewb.[{£/2 + x T )y/sjD] 



x sinh[(^/2 - x r $) y/s/D] , x T < x T ,o 



x < 



(E4) 



sinh[£y/s/D] 
-f LX [% S mh[{£/2 



x T )yf7/D] 
x sinh[(£/2 + XT,o)y/s/D], xt > xt,o 



1 



ssmh[£y/s/D] 
smh[£y/s/D] 
-fmf§ sinh[(£/2 • 



Xr 



x sinh[(^/2 — xt, o)y/ s/D], xt < ir,o 



x < 



(E5) 



(1 - {xt - XT,o)f R ) smh[£y/I/D] 
sinh[(£/2 - zr)^7^] 
x sinh[(^/2 + XT,o)y/s/D], xt>xt,q- 



The remaining three integrals follow from the normaliza- 
tion conditions (arguments left implicit): 



1 



n+n 



(E6) 



Laplace inversion of the above relationships, using e.g. 
residue calculus [43|], gives Eq. (f2T))) . In the following 
subsections we give asymptotic results in the (1) long 
and (2) short time limit for the expressions above. 



1. Long-time behavior 

The long-time behavior of Eqs. (|E1[) - (|E16[) is obtained 
from a series expansion for £y/s/D <C 1, and reads (ar- 
guments left implicit) 



1 



s \2 



(E7) 



The inverse transforms are found from C 1 (s 1 ) = 1 44 1. 



Short-time behavior 



Short times is defined here as {£ ± xt)\Js/D, (£ ± 
xt ,0) V s I D 3> 1, i.e. times shorter than the time it takes 
to diffuse across the entire box. The short time behavior 
of Eqs. (|ET|) - (|ET6)) is given by 



2s e 

2 C 



's/D 

■XT,o)y/s/D 



XT < XT,0 
XT > X T ,0 

(E8) 
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-{xt,o—xt}\/ s/D 



ft r -(xT-x r .a) 



2s ' 



s/D 



XT < XT,0 
XT > XT,0 

(E9) 




1 fl - M /R 

s \ 2 V « 

\{l-{xT-X T , a )f, 



(ir,o-ir)V s / D 



V X T < X T fl 
XT > XT,0 

(E12) 

The remaining integrals follow from Eq. (|E16|) . 

Equations (|E9|) - (|ET2|) can be inverted exactly into time 
domain. Using standard formulas (44j and introducing 



Sr D_ p -( x i — xt,o)\/s/D 
2 



XT ~ XT,0 



(E13) 



leads to 



^(l-erftj), V*(*) = y(l + erftj) 
i(l+erf^), ^(t) = i(l-erf» 7 ) 



^(t) = 1 - 



f R 4Dt 



e v + \/7rrj(erf?7 — 1) 
e _r ' 2 + v/Tr^Cerfr? + 1) 



(E14) 



where erfz is the error- function [44j |. The expressions 
above are valid for times such that 4Dt/(£/2±XT,o) 2 1 
and ADt/(l/2 ± ir) 2 < 1. Expanding the result in Eq. 
(|E14|) for small 77 and using the normalization conditions 
Eq. (|El6|) . yields 



^(*) 



1 - 



1 - Vttt]) , 



/4Dt , r- \ 



/fl. 41)/ 



and (arguments left implicit) 



(E15) 



(E16) 



The limit where i — > limit (s — > 00) is most conveniently 
found from Eqs. (|E1[) - (|E16|) which, in combination with 
£ _1 (s _1 ) = 1, reads 



^z(t->0) = 



« -> 0) = 



0, XT > XT.Q 



0, sr < £r,o 
/r, st > xt,o 



(E17) 



(E18) 



(E19) 



^(i - 0) = ^ } - {XT >° ~ Xr)/L ' ^ | xr '° (E20) 

1, ^ X'J' Q 



^«(*-0) = (J' , w ^sr.o (E21) 

v flV 7 \ 1 - far - X T fi)fR, x T > x T , a v 7 

Again, the remaining integrals follow from the normal- 
ization condition (|E16[) . 



APPENDIX F: MACROSCOPIC DYNAMICS - 
THE DYNAMIC STRUCTURE FACTOR AND 
CENTER-OF-MASS MOTION 

Macroscopic quantities for a single-file system is the 
same as for a system consisting of non-interacting par- 
ticles [l9|, [I?]] • By "macroscopic" we here refer to any 
quantity which is invariant under Xi — > Xj, i.e. any one 
which do not " notice" if two particle are interchanged in 
the system. In this Appendix we explicitly evaluate two 
such macroscopic quantities: (1) the dynamic structure 
factor and (2) the PDF for the center-of-mass coordinate. 
We demonstrate that indeed they agree with known re- 
sults for non- interacting systems. 



1. Dynamic structure factor 

The dynamic structure factor is tractable via scat- 
tering experiments (see e.g. Ref. (Sf|), and is widely 
used in condensed matter physics and crystallography. 
Considering a set of noise-driven stochastic trajectories 
Xi(t), ...,Xjy(t), the dynamic structure factor is (6?| 



iQ(Xi(t)-X jl0 ) 



(Fl) 



where the brackets denote an average over different re- 
alizations of the noise. In terms of the jVPDF and the 
equilibrium iVPDF (see Sec. HI}, it is given by 

S(Q,t) = i / dxi---dxN [ dxifi ■ ■ ■ dxN,o 
iV Jtz Jk 

x J2 etQiX '~ X3,o) ' P (x> t \xo)V cq (x ) (F2) 

i,3 
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where we averaged over the (equilibrium) initial posi- 
tions. Since the integrand above is invariant under 
Xi «-> Xj we can apply the technique explained in Ap- 
pendix [C] to extend the integrations over coordinates as 
well as initial positions to [—£/2,£/2], which leads to 



where S(z) = (2n) 1 da e laz was used. Inserting the 
results above into Eq. (|F6j) gives 



S(Q,t) = e 



- p -DQ 2 t 



(F9) 



S(Q,t) 



N Nl £ N 1 ' 



1/2 



dx 2 

-1/2 J -1/2 
1/2 j-l/2 

dxifi j dx 2 .o ■ ■ ■ 

-1/2 J-l/2 



dx 



N 



dx N fi 



(F3) 



which is in agreement with standard results for non- 
interacting particles in one dimension (see e.g. Ref. 68]). 



2. Center-of-mass dynamics 

The probability density function for the center-of-mass 
(CM) coordinate X is given by 



where also Eq. (j4j was used. Inserting explicitly the 
iVPDF from Eq. (16]) yields 



P(X,t\x ) 



dxidx2 ■ ■ ■ dxwS ( X — 



n 



X\ + x 2 



■ x N 



N 



i i i r'* 
S (Q>*) = NmWj_ e/2 dxi 

/2 

dxin 



dx 



2 ■ ■ ■ 



dx 



N 



xP(x,t\x ), 



(F10) 



dx 



2,0 ■ 



dxjsj ( 



N 



i=l 



xip(xi,xifi] t)ip(x 2 ,x 2fl ; t) ■ ■ ■ ip(x N , x Nj0 ; t) 



where the 7VPDF P(x, t\xo) is given in Eq. (fl2)) . and the 
integration region 1Z is defined in Eq. ([3]) . Since the inte- 
grand is invariant under it is possible to extend 
the integration region to Xj £ [—1/2,1/2] (provided we 
divide by Nl), as was done in previous subsection. Also, 
using the integral representation of the (5-function from 
previous subsection, Eq. (|F10[) can be rewritten as 



remaining AH — 1 terms 



(F4) 



Defining 
C(Q) = 



P(X,t\x ) 



x exp 



dx, 



1 

m 
x - 



dQ 
2ir 



dx\ ■ 

oo *"■ J-l/2 
X\ + X 2 + ... + Xn 

N 



dx 



N 



P(x, t\x*o). 
(Fll) 



Combining this equation with Eqs. (fT2|) and (fl3|) . and 



= p J_ t 2 dXl J 2 dx 3'° J dxk J 2 dxi <° defining g{k ,Q) = ff /2 dxi ^^ 3 -Q/N) gives 



x^( Xi , x jt0 ;t)il>(x k ,xt, ; t)e^ x '~ Xl - o) (F5) 

and noticing that all the Nl terms corresponding to per- 
mutations of the initial particle positions in Eq. (|F4|) 
give the same contribution, leads to 



P(X,t\x ) 





r°° dk! 


r°° dk N 


Lao 27T 


Loo 2tt " 


J-oo 27T 



xe- D & + ~ +i & t 4>(ki,xi, ) ■ ■■4>{k N ,x N , ) 
xg(h,Q)---g(k N ,Q). (F12) 



S(Q,t) = C(Q)C(0) 



N-l 



(N - 1)F(Q)C(0) . (F6) Integrating over jb leads to 



This result is identical to that of non- interacting particles 
(obtained e.g. by putting, SV, = 0, in the Bethe Ansatz 
Eq. P^|) and allowing the particles to be in the full phase 
space xj £ [— 1/2, £/2]). The expression above is simpli- 
fied by noticing that C(0) = 1, due to the normalization 
of ip(xi,Xjfi;t). 

As an example, we consider the case i — > oo and 
A = where ip(xi,Xjfi-,t) = {AirDt)^ 1 ^ 2 exp[— (xj — 
x j,o) 2 /(4Di)] can be used. A straightforward calculation 
of Eq. (IF5p for this case gives 

C(Q) = e-°Q 2t 
F{Q) = ^S(Q)e~ D ^, 



P(X, t\x ) = I ^-e lQX h{Q, xi, ) • • • h(Q, x Nj0 ), 
, 2-ir 



where 



(F13) 



f°° (ft., 2 

h(Q,xjo) = / tt Z -e~ Dfc 'V(*j,*jo)ff(*j.Q) 

J-oo 27T 

" e/2 f°° dk 

dxi I ^l e -Dk^t 

2tt 



(F7) 
(F8) 



x^(fc„x,,o)e Ml(fc ^ Q/Af) 



(F14) 



The result in Eq. (|F13[) is identical to that of non- 
interacting particles, as can be shown straightforwardly 
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by setting Sij = in Eq. (fT2|) and not restricting the 
particles to the phase space region 1Z. 

As an example, Eq. (|F13|) is evaluated for an infinite 
system (£ — > oo) and A = where <j)(kj, £j o) = e~ lkjXj -°. 
For this case Eq. (|F14[) becomes 



h{Q,Xjfl) = exp 



Q 2 Dt iXjfiQ 



(F15) 



N 2 N 

which when inserted in Eq. (|F13[) and integrated over Q 



gives the well known Gaussian for the CM 



P(X,t\x ) 



1 



(AnDcMt) 1 / 2 



exp 



yO \2 
^CMJ 



4DcMt 



(F16) 

where Xq M — fj^i^i x i,o) /N and D CM — D/N denote 
the CM initial position, and diffusion constant, respec- 
tively. 
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